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A method of computer simulation of Time-Resolved X-ray Diffraction (TRXD) in asymmetric 
Laue (transmission) geometry with an arbitrary propagating strain perpendicular to the crystal 
surface is presented. We present two case studies for possible strain generation by short-pulse 
laser irradiation: (i) a thermoelastic-like analytic mode^; (ii) a numerical model including effects 
of electron-hole diffusion, Auger recombination, deformation potential and thermal diffusion. A 
comparison with recent experimental results- is also presented. 

PACS numbers: 07.05.Tp, 61.10.-i, 78.47.+p 
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I. INTRODUCTION 

When a semiconductor crystal is irradiated by photons 
of above band-gap energy, electrons are excited from the 
valence band to the conduction band. If the light source 
is a subpicosecond laser, the electron-hole plasma density 
can easily reach 10 20 cm -3 — enough for the lattice spac- 
ing changes due to the deformation potential to be signif- 
icant and for the fast diffusion of the electron-hole plasma 
into the crystal bulk to be appreciable. The time-scale of 
the fastest energy transfer processes from the electrons 
to the lattice is of the order of picoseconds or less, which 
is significantly shorter than the hydrodynamic response 
time of the crystal. The resultant stress is therefore re- 
lieved by surface expansion and, by Newton's third law, 
a bipolar compression wave propagating into the crystal. 
These effects have been studied by optical methods 3-4 - 5 - 6 
and, more recently, as test-cases for the field of Time- 
Resolved X-ray Diffraction (TRXD) in Bragg (reflection) 
geometry li&^iiSiiiiiSii^ However, although x-ray diffrac- 
tion in Bragg geometry allows the direct study of the 
structural changes to the lattice — as opposed to the elec- 
tronic effects in optical studies — the very nature of Bragg 
geometry x-ray diffraction with perfect crystals only al- 
lows the study of the surface region of the crystal. By 
utilising the anomalous absorption effect in Laue (trans- 
mission) geometry, recent experiments 2 ^ have enabled 
the probing of structural changes throughout the bulk 
of the crystal on a sub-nanosecond time-scale, giving a 
much clearer picture of the phenomena under considera- 
tion. 

In Bragg geometry it is only possible to study the 
structural changes of crystalline matter in a layer with 
thickness comparable to the shorter of the x-ray extinc- 



tion depth (~ (im) for strong diffraction or the absorp- 
tion depth for weak diffraction. 9 In Laue geometry it is 
possible to study the whole bulk of the crystal due to 
anomalous absorption, known as the Borrmann effect. 
In dynamical diffraction theory, the x-ray electromag- 
netic field inside the crystal is resolved into two inde- 
pendent eigensolutions — one with nodes on the lattice 
planes, and the other with anti-nodes. This creates one 
solution with reduced absorption, and another with en- 
hanced absorption. These are known respectively as the 
a and f3 branches. It is the a branch solution that can 
propagate for many extinction depths without significant 
loss — probing the entire depth of the crystal^ A distur- 
bance in the crystal lattice causes energy to be trans- 
ferred between the two branches, an effect known as inter- 
branch scattering* 1 ^ When the x-ray beam exits the crys- 
tal, it is again expressed as the free-space solutions: the 
(forward-diffracted) and h (deflected-diffracted) beams, 
which are linear combinations of the a and f3 branches. 
If the ratio and/or relative phases of the branches has 
changed, the energy partitioning between the beams will 
be affected, often to a significant degreed 



Using this method it is possible to study the mecha- 
nism of energy transfer when a single crystal target is ir- 
radiated by a femtosecond laser pulse. We will show that 
it is possible to distinguish between a model which as- 
sumes an instantaneous energy transfer into the lattice, 
and one which models thermal diffusion, electron-hole 
plasma diffusion, the effect of the deformation potential, 
and energy transfer to the lattice via Auger recombina- 
tion. 
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FIG. 1: Layer approximation. (1) At the boundaries between 
layers, the a and f3 branches are calculated (Eg. 111(1 . (2) The 
beam is then propagated through the layer (Eq. |1J . 



II. THEORY 

The fundamental assumption in this method of diffrac- 
tion simulation is that the strain field in the crystal is 
one-dimensional, parallel to the surface normal, and can 
be approximated by a constant strain in many laminas, 
each one parallel to the surface of the crystal. The de- 
gree to which this represents the experimental situation 
is discussed in mill The electric displacement field ampli- 
tude is then calculated as it passes through these laminas, 
taking into account the effect of the strain in each. 



These vector equations relate the physically realisable 
values of T)o,h and can be solved for the incident and 
diffracted beam amplitude as a function of depth into 
a crystal layer giving equations of the form (for a- 
polarization) 
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are the wave- vector shifts of the a and (3 branches (nor- 
mal to the crystal surface — see Fig. H), D a and Dp are 
the amplitudes of the a and (3 branches, 
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A. Dynamical Diffraction Theory 

By solving the wave equation 

curl curl [(1 - x)D(r)] = 47r 2 fc 2 D(r) (1) 
in a medium with a periodic susceptibility 
X (r)=^x^- Mh - r , 

h 

using a displacement field 

D(r) =^D ft e^ ik -- r , 

h 

where k/j = ko + h is the wave- vector of the h-ih. compo- 
nent inside the crystal and both are summed over fourier 
components in the reciprocal lattice vectors h, we obtain 
the fundamental equations of dynamical theoryii 

J2 Xh-h> [(k& • Dv)k& - k 2 h T> h ,] = (k 2 - kl)T> h . (2) 

h' 

These can be further simplified by the assumption that 
there are only two strong x-ray wave-fields (incident and 
diffracted: D(r) = D e~ 2Triko ' r + D ft e- 2Trik "- r ), leading 
to the two-beam dispersion relation 

X - h [(k • D,0k o - klD h ] = [k 2 - k 2 (l - xo)] D , (3a) 
Xh [(k,, • D )k h - fc 2 D ] = [k 2 - k 2 (l - xo)] D ft . (3b) 
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A9 = 9 



(9) 



where 70, "fh are the direction cosines of the incoming and 
outgoing beams respectively, in relation to the surface 
normal of the crystal n; £ aj( 3 are the amplitude ratios 
Dh/Do; 9 is the actual glancing angle of incidence relative 
to the diffracting planes; and 9b is the Bragg glancing 
angled One can also define a dimensionless deviation 
parameter 77, defined as 



A9 - A6 n 



(10) 



where A9 os is the refraction shift of the diffraction peak 
and 26 os is the Darwin widths This can be shown to 
be equal to — -^= in terms of the quantities shown above. 
The deviation parameter can also be expressed as a func- 
tion of wavelength separation from the Bragg wavelegth 
A A = A — A B , using the relation AA = (Ab cot 9b) A9 
derived from Bragg's law. This is the form that will be 
used for analysing the experimental data. However, it is 
easier to visualise the effects of strain (see the following 
section) when it is considered as a function of angle. 

Eq.^can be solved at the boundary between two layers 
(n • r = 0) to find D a and Dp as functions of Do and Dh 
in the previous layer: 
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where the superscripts indicate the layer in question. It 
should be noted that k a _p and £ Qj/ 3 are all functions of 
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FIG. 2: Strained crystals. Strain (ej_ = — ) alters the rotation 
(S(p) and the separation (Sd) of lattice planes in the crystal. 
Note that there are two possible orientations for the x-ray 
beams. 



77, and it is through this parameter that the strain is 
incorporated. Therefore, in Eq. 1111 care must be taken 
that the £ Qi( 3 are evaluated in the correct layer. We may 
do so by considering each layer as a separate crystal, 
with the beams D^ h exiting crystal (i) and then entering 
crystal (i + 1) without change. It is then easy to see that 
the correct layer in which to evaluate these quantities is 
layer (i + 1). 



The change in the Bragg angle can be calculated using 
the differential form of Bragg's law, 

5d 

59b = — r tan 
a 

combined with Eq. 1131 giving 

-59b = £j_ sin 2 </>tan6> B . (15) 

We can now include both of these effects in the angu- 
lar deviation from the Bragg angle to give the 'effective 
misorientation' 

(9-9 B ) -> (6>-6> B ) + (sin 2 0tan6lBTsin0cos0)e ± . (16) 

This can be used in every lamina, with the appropri- 
ate strain, in order to calculate the time-resolved rock- 
ing curves from the strained crystal. In order for this 
lamellar approximation can be considered accurate, each 
lamina must be sufficiently thin that the change in the 
effective misorientation is small within it. In practice this 
can be achieved by using a variable layer thickness, with 
a maximum allowed strain change across a layer. The 
maximum strain change within a layer can then be re- 
duced until the results obtained remain constant. This 
was achieved with a maximum strain change of 1 x 10~ 8 
within any layer, for the geometry under consideration 
in this article. 



III. STRAIN MODELS 



B. Incorporating Strain 

Strain is incorporated into the model through the de- 
viation parameter 77 (or more specifically, its real part 
i] r ), which can be expressed as a function of the angu- 
lar deviation from the Bragg angle A9 (Eq. I1U|) . Strain 
affects this angular deviation in two ways: (i) rotating 
the planes (changing 9) and (ii) altering the separation 
(changing 9 B )- The method presented here is similar to 
that of earlier work in Bragg geometry ji&iiLSl but is now 
presented for the Laue geometry case. 

From Fig. [3 we can formulate the following relations 
by taking differentials of the relationships between d, x, 
y and 4>. 

5(f) = — ej_ sin 4> cos <j>, (12) 
5d o 

— = e_Lsin <p, (13) 

where e± = & is the strain perpendicular to the surface. 

The change in diffraction angle 59 is equal to ±5(f) de- 
pending on the orientation of the incoming and outgoing 
beams (6\ or 9 2 in Fig. |21 for positive and negative <f> re- 
spectively). The change in diffraction angle due to the 
planar rotation is therefore 

56 = =Fe_i_ sin <j> cos <j>. (14) 



As an example of this technique, we use the case of 
a Germanium single crystal irradiated by a femtosec- 
ond pulse of near infra-red laser radiation. Experi- 
ments using Laue geometry time-resolved x-ray diffrac- 
tion on such a sample have recently been carried out. 2 
This case has also been studied experimentally using 
ultrafast reflectivitji 3 ^^ and x-ray Bragg scattering 
techniques, 7 ! 8 ! 9 ! 10 ! 11 . 12 ! 13 

In this experimental set-up, the short laser pulse 
quickly heats the surface at constant volume, generat- 
ing a thermal stress. This causes the surface to expand 
and, by Newton's third law, launches an acoustic pulse 
into the crystal. In both cases considered here, the lat- 
eral size of the laser spot is assumed to be much greater 
than the laser absorption depth. As a result, the strain 
generated can be assumed to be one dimensional^ i.e. 
the atomic displacement only varying as a function of 
depth into the crystal. 

We will be comparing two models of the strain pro- 
duced in this system. The first model takes into account 
the time-scales of the processes by which the energy is 
transferred to the lattice. The second is a simplified case 
with an analytic solution, introduced by Thomsen et alX 
It assumes instantaneous transfer of energy from the laser 
into the lattice and no diffusion. 

When a short laser pulse is incident on the crystal, the 
following energy transfer processes occur (see Fig. |3J): 
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FIG. 3: Band structure of Germanium. The transfer of en- 
ergy to the lattice involves: (1) Laser excitation. (2) Fast 
relaxation to band edges. (3) Auger recombination. 



1 . The laser excites electrons from the valence band to 
the conduction band of the semiconductor, creating 
an electron-hole plasma. 

2. The electrons and holes quickly relax to the band 
edges, transferring their energy to the lattice. This 
energy transfer is assumed to be instantaneous on 
the time-scale considered (< lps). 

3. The electrons and holes then recombine (by Auger 
recombination — an e-e-h and e-h-h process), trans- 
ferring the rest of the energy to the latticed This 
takes significantly longer (typically 1 ns to 1 /us) , al- 
lowing the carriers to diffuse appreciably into the 
crystal before recombining. 

The electron-hole plasma and the thermal phonons in 
the lattice both diffuse into the crystal bulk, with sep- 
arate diffusion constants. They both cause a change in 
the lattice spacing, through the deformation potential 
and thermal expansion respectively. 

The laser energy is deposited in certain area of the 
crystal with a 1/e absorption depth £ and an absorbed 
fluence Q. This gives initial conditions of the electron- 
hole plasma and the lattice temperature as 



n(z,t = 0) 



(17) 



T(z,t = 0) = ^— %Ln(z,t = 0), (18) 

where n is the electron-hole plasma density, E p is the 
energy of the laser photons, E g is the indirect band gap 
and Ci is the lattice heat capacity per unit volume. 

The electron-hole plasma obeys a diffusion equation 
with a sink term for Auger recombination 



dn d 2 n 

~dt = p d^ 



(19) 



where D p is the plasma diffusion constant and A is the 
Auger recombination rate. The energy from the Auger 



recombination is then transferred to the lattice, which 
obeys a diffusion equation with a corresponding source 
term 



dt * dz 2 



where Dt is the thermal diffusion constant. 
The equilibrium strain is then 

e e (z, t) = a t T(z, t) + a p n(z, t), 



(20) 



(21) 



where at is the thermal expansivity (= <E>/3, $ is a factor 
to take into account the ID nature of the strain, (3 is 
the linear expansion coefficient) and a p is an electronic 
contribution to the strain associated with the deforma- 
tion potential (= 9(log a)/dn, a is the equilibrium lattice 
constant). Any change in the equilibrium strain produces 
forward and backward propagating waves, which can be 
calculated by integrating Eq. |^as follows: 

e+(z,t) = ~J ^- t [e e (z-vt',t-t')]dt', (22a) 



e-(z,t) 



\j §1 M* + <*-*')]< (22b) 



where v is speed of longitudinal sound in the crystal. The 
free-surface boundary condition of the stress, 033 (z = 
0, t) = 0, is ensured by defining 

e e (-z, t) = -e e (z, t), z > 0. 

The total strain is then the sum of the equilibrium and 
the forward and backward going strain waves. 



e(z,t) = e e (z,t) +e+(z,t) +e-(z,t). 



(23) 



The Thomsen strain mode! 1 in its most simple form 
assumes instantaneous transfer of energy into the lat- 
tice and no diffusion. In the time-scale of interest only 
a fraction (E p — E g )/E p of the absorbed laser energy 
is transferred to the lattice (the rest remaining in the 
electron- hole plasma). This gives a total strain of 



(CiE p 
1 r 



(*+«*)/C +e -|»-trt|/i sgn(z _ ut) 



(24) 



Note that this is the sum of a non-evolving exponential 
surface strain, a forward going wave T(z — vt) and a 
backward going wave Q(z + vt). 

The numerical model produces an identical strain to 
the above analytic formula, if there is no deformation 
potential or Auger recombination (a p = and A = 0). 
A comparison of the strains produced by the numerical 
and Thomsen strain models is shown in Fig. 01 
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FIG. 4: Strain profile evolution — comparison of numerical and 
Thomsen models. This figure shows a time sequence of strain 
snapshots taken at ~ 40 ps intervals (v = 4915 ms -1 ) for an 
absorbed fluence of Q — 4mJ/cm 2 . The physical parame- 
ters used are the literature values given in JVTl The strain 
in the numerical model comprises a diffusing surface compo- 
nent and an asymmetric travelling wave. The asymmetry is 
caused by the e-h diffusion, the higher peak strain by the 
deformation potential, and the surface decay by a combina- 
tion of thermal and e-h diffusion. The strain in the Thomsen 
model is made up of an exponential surface component and 
an anti-symmetric travelling wave. The travelling portions of 
the strain move into the crystal at the speed of sound. 



IV. SIMULATION 
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FIG. 5: Simulated time-resolved rocking curves. The and 
h beam rocking curves plotted as a function of time for the 
numerical strain model with an absorbed pump fluence of 
Q = 17mJ/cm 2 . 
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FIG. 6: Experimenal setup. The h beam is collected by an 
avalanche photo-diode (APD) giving a time resolution of the 
synchrotron bunch length, the beam is resolved by a streak 
camera giving picosecond resolution. 



The diffraction simulation works as follows: Over a 
number of time-steps, the strain profile (analytic or nu- 
merical) is calculated. For the numerical method, the 
diffusion equations (|19l and l2"U)l are solved using a Crank- 
Nicholson scheme, and the integrals in Eq. [22 are cal- 
culated by finite differencing then summing. The x-rays 
are then propagated through the crystal over the rj r range 
of interest. This is done using the propagation matrices 
f flll A|) . After the final layer, the amplitudes for the 
and h beams are multiplied by their complex conjugates, 
giving the rocking curves (Fig. [SJ) . The resulting time- 
resolved rocking curves are then integrated over r\ r and 
normalised such that the value for an unstrained crystal 
is unity, for comparison with experimental results. 

It is also possible to output the beam intensities 
through the bulk of the crystal: either as the and h 
beams, showing Pendellosung oscillations 15 ; or as the 
a and f3 branches, showing the transfer of energy be- 
tween the two as the strain moves through the crystal 
(see Fig.EJ. 

V. EXPERIMENT 

The experimental data was taken at the 7-ID undula- 
tor beamline at the Advanced Photon Source. The x-ray 
beam energy was lOkeV with a 1.4 x 1CP 4 fractional en- 



ergy spread (larger than the rocking curve width of the 
crystal) and negligible beam divergence. The sample was 
a 280 /im thick (001) Germanium single crystal, oriented 
to diffract from the 202 planes (Fig.©. The x-ray beam 
was masked by tantalum slits giving a beam size on the 
crystal of 400 /im x 400 /im. The strain pulse was pro- 
duced by exciting the output face of the crystal with 
sub-lOOfs, 800 nm laser pulses. The diffraction simula- 
tion was set up to match these experimental conditions. 

VI. RESULTS AND DISCUSSION 

Data taken experimentally 2 is shown in Fig.[7{a). Os- 
cillations are seen in the h beam amplitude, with the time 
evolution of these varying with pump fluence. The phase 
of the oscillations is calculated by fitting curves of the 
form 

I(t) = ^a,e~^ sin(wii - fc), (25) 

i 

starting at t = 1 ns (to ignore the initial rise) . The size 
of the initial rise or fall is measured by taking the nor- 
malised intensity at t — 200 ps. The phase <f> and the ini- 
tial rise or fall (\D h (t = 200ps)| 2 / \D h (t = 0ps)| 2 ) of the 
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FIG. 7: Time-resolved integrated reflectivity. Plots are for diffracted beam (|D;,| 2 ). (a) Experimental results, (b) Thomsen 
model with £ = 200 nm. (c) Thomsen model with £ = 1 /im. (d) Numerical model, Dt = 0.35cm 2 /s, D p — 65cm 2 /s, 
A = 1.1 x 10~ 31 cm 6 /s, at = 10" 5 K _1 , a p = 1.3 x 10~ 24 cm 3 , £ = 200 nm, E g = 0.67 eV, C ; = 1.7J/Kcm 3 and v = 4915ms" 1 . 
(e) Numerical model with no deformation potential, as (d) but with a p = 0. (f ) and (g) Numerical model with increased Auger 
recombination, as (d) but with A — 5.5 x 10 -31 cm 6 /s and 11 x 10 -31 cm 6 /s respectively. Note that modelling assumes that 
75% of the incident laser energy is absorbed. 



measured and simulated reflectivity is shown in Fig. [S] 
The curves for measured reflectivity assume 75% of the 
laser energy is absorbed by the crystal. 

Simulations were run with the Thomsen strain model, 
using literature values for physical parameters (C = 
200 nm, E„ = 0.67 eV, <f>/3 = 1 X 10 _5 K _1 , Q = 
1.7 J/Kcm 3 ). The results are shown in Figs. [7| and Efb). 
These match well with the phase change of the oscilla- 
tions (to approx 10%), but show a very different overall 



form and don't match the initial rise. If the laser absorp- 
tion depth £ is increased to 1/xm (Figs. Handle)), the 
overall form of the higher fluence curves shows a better 
match. However, the phase of the oscillations doesn't 
show the expected behaviour. 

Using the more detailed numerical model, again with 
literature values for physical parameters (C = 200 nm, 
a t = 1 x lO^KT 1 , a p = 1.3 x 10" 24 cm 3 , A = 1.1 x 
10" 31 cm 6 /s, D t = 0.35cm 2 /s, D p = 65cm 2 /s), the 
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FIG. 8: Effects of changing the deformation potential and the 
plasma diffusion constant on (i) the phase of the oscillations 
and (ii) the initial rise or fall of |-Dh| 2 . The phase is fitted to 
curves of the form shown in Eq. 1251 with the data starting at 
t — 1 ns to ignore the initial rise. Fitted to: (a) Experimental 
data (assuming 75% absorption), (b) Thomsen strain model, 
£ = 200 nm. (c) Thomsen model, £ = 1 /im. (d) Numerical 
strain model, a t = 1 x 10~ 5 K" 1 , a p = 1.3 x 10~ 24 cm 3 , A = 
1.1 x 10 -31 cm 6 /s. (e) Numerical model with no deformation 
potential, as (d) but with a p = 0. (f) and (g) Numerical 
model with increased Auger recombination, as (d) but with 
A = 5.5 x 10~ 31 cm 6 /s and 11 x 10~ 31 cm 6 /s respectively. 



curves in Figs. [7| and |HJd) were calculated. These match 
the fluence dependence of the experiment in the over- 
all form of the curves and the initial rise. However, the 
phase of the oscillations only decreases at about two- 
thirds the rate of the experimentally measured phase. 
At the highest fluences the oscillations in D) x are of such 
small magnitude that it was impossible to fit the curves 
numerically. 

If the deformation potential is taken out of the nu- 
merical model (by setting a p — 0), the correct fluence 
dependence of the phase is obtained (to approximately 
10%), but the behaviour at early times (ie. the initial 
rise or fall) is no longer well matched (Figs. and |SJe)). 
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FIG. 9: The //-integrated interior solutions \D a \ (dot-dashed 
line) and \Dp\ 2 (dashed line) as a function of depth into the 
crystal at 1.25 ns after laser excitation, for an absorbed pump 
fluence of Q = 17mJ/cm 2 . Also shown is the strain distri- 
bution (solid line). Strain front is moving from right to left. 
X-ray beams exit crystal on right. 



These last two observations would imply that there is 
a direct electronic contribution to the strain (through 
the deformation potential) at earlier times, which is 
greatly reduced at later times. One possible mecha- 
nism for this is that the Auger recombination rate is 
larger than expected. Simulations were run with the re- 
combination rate increased five-fold and ten-fold (A = 
5.5 x 10~ 31 cm 6 /s and 11 x 10 _31 cm 6 /s respectively). 
The resulting curves (Figs. [7| and 03 (f) and (g) re- 
spectively) show good agreement with the experimental 
phase, and a reasonable agreement with early time be- 
haviour. 

We were unable to find an excellent match with any 
combination of parameters, but this particular set gave 
the best of those tried. The exact mechanism for a larger 
than expected Auger recombination rate is unknown. 
One possible explanation is that the higher temperature 
of the lattice at early times activates other possible re- 
combination pathways (leading to a temperature depen- 
dence to the recombination rate). It is thought that the 
main reason for the recombination rate being so low in 
Germanium is that it requires phonon activation^ So, 
at the higher temperatures existing at early times, this 
could cause faster than expected recombination. More 
detailed modelling of the strain, taking this possible non- 
linearity into account, could lead to better reproduction 
of the experimental data. 

In an earlier paper ^ the authors posited that the basic 
physics behind the oscillations visible in the reintegrated 
rocking curves (Fig. 0) was due energy being transferred 




FIG. 10: 2D Region (drawn to scale). The incoming Do beam 
has a width of 400 /im, the crystal is 280 /im thick. In the '2D 
region' shown, the exit surface of the crystal is not influenced 
by the whole incoming beam. 



from the a branch to the strongly absorbed (3 branch at 
a disturbance in the lattice. Using the diffraction sim- 
ulation it is possible to look at how the strain transfers 
this energy between the solutions. Fig. shows that en- 
ergy is transferred from the a branch to the (3 branch 
when the strain gradient is negative, and is then quickly 
transferred back at the positive strain discontinuity. This 
energy transfer effect is reversed for the opposite asym- 
metry (202). 

It might initially be expected that analysis of these ex- 
perimental results would need a 2D diffraction model. 
This is because, as shown in Fig. I1UI for this experi- 
ment only a small fraction of the exit surface is influ- 
enced by the entire incoming x-ray beam — the area of 
influence being the inverted Borrmann triangle from the 
exit points The 2D Takagi-Taupin equations normally 
reduce to a case only dependent on the depth into the 
crystal for a one dimensional strain, as in this case^ 
However, this also has the implicit assumption that the 
lateral extent of the x-ray beam is large enough for its 
size to be unimportant. 2D simulations were run, us- 
ing the 'half-step derivative' numerical solution^ Sam- 
ple results are shown in Fig. 1111 For comparison with 
the ID results, the intensity distribution was summed 
over the exit surface of the crystal, and then over the 
rocking curve. The time-resolved integrated reflectivi- 
ties obtained only differed by 0.25% (RMS) from the ID 
diffraction model. The reason for this is that close to 
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o 1 
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FIG. 11: 2D simulation results (drawn to scale). [D^f within 
the crystal for (a) rj r = —1, (b) r\ r = 0, (c) r\ r = +1. 



the rocking curve peak, the D^ ^ beams are not 'eigen- 
solutions' inside the crystal — the D a ^ solutions are, due 
to the strong coupling between the beams. We would 
therefore only expect differences to become apparent in 
the rocking curve wings, where the coupling is less strong. 
In the experimental conditions, where the data obtained 
is the Ty-integrated rocking curve, the great majority of 
the signal comes from the rocking curve peak, completely 
swamping any small changes in the wings. 



VII. CONCLUSION 

We have presented a method for the calculation of 
time-resolved rocking curves for X-ray diffraction in Laue 
(transmission) geometry. Such diffraction studies are an 
important technique for the study of coherent strain in 
crystals beyond the extinction depth. Even without the 
benefit of rocking curve resolution of the diffracted x- 
ray beams, it is possible to obtain valuable information 
about the form of the strain profile in the crystal, and to 
observe the effects of the mechanisms of the ultrafast en- 
ergy transfer processes at work. As there is no method to 
analytically compute the strain in a crystal from a time- 
resolved rocking curve, numerical simulations are key for 
the understanding of experimental results obtained. 
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